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Recently, much attention has been devoted to the study of two-dimensional walk 
models and their relations with the one-dimensional driven- diffusive systems [HOE]. 
The reason is that the partition function of some of these walk models can be related to 
those of some of the one-dimensional out-of-equilibrium driven- diffusive systems with 
open boundaries. One of the most simple and important system of this type is the 
Asymmetric Simple Exclusion Process (ASEP) with open boundaries [I]. 
The ASEP is a widely studied driven- diffusive system which reveals a wealth of 
interesting critical phenomena. The ASEP is defined on a one- dimensional discrete 
lattice. The classical particles are allowed to enter the lattice only from the left 
boundary and leave the lattice only from the right boundary. The particles also hop 
unidirectionally toward the right boundary. It is known that a single product shock 
measure has a simple random walk dynamics in this system, without any constraints 
on the parameters of the system, provided that we consider a discrete-time updating 
scheme. In this case the particles are injected into the leftmost lattice site with the 
probability a if the target site is empty. They are also extracted from the rightmost 
lattice site with the probability (5 if it is already occupied. We assume that in the 
bulk of the lattice a particle at a given site k deterministically move to the site k + 1 
provided that it is not occupied by other particles, which means that double occupation 
is prohibited. 

In this paper we consider the ASEP with open boundaries in the sublattice-parallel 
update [5], E]. We are looking for a two-dimensional walk model which shares similar 
critical properties with the ASEP. It turns out that there exists such a two-dimensional 
walk model as we will describe it later. These systems have equal partition functions and 
the physical quantities in these systems are related through a simple transformation. 
The ASEP is defined on a lattice of length 2n. The time evolution of the probability 
distribution vector \P(t)) is given by the following master equation: 

T\P(t)) = \P(t + l)). (1) 

The sublattice-parallel updating scheme is defined as follows: We divide the bulk 
dynamics consists into two half time steps. In the first half time step the pairs of 
neighboring sites (2k, 2k + 1) for k — 1, • • • , n — 1 and also the first and the last lattice 
sites are updated. In the second half time step the pairs of neighboring sites (2k — 1, 2k) 
for k = 1, ■ ■ ■ ,n are updated. Considering this definition, the transfer matrix T in (TjQ) 
is given by the multiplication of two factors T = T 2 7i defined as: 

t x = c®r®...®r®n = c®r® {n - l) ®n 

T 2 = T®T®...®T = T® n 
where T, C and 1Z are given by: 
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The matrix T is written in the basis (00, 01, 10, 11) when stands for an empty lattice 
site and 1 stands for an occupied lattice site. The matrices £ and 1Z are also written in 
the basis (0, 1). 

Following [8] we define two different types of product shock measures |/x 2fc ) (k = 1, ■ ■ • , n) 
and |/i2fc+i) (k — 0, • • • , n) as follows: 





V Pi J \ Vi J \ 

2/c+l 

Note that the position of the shock in \[i2k) lies between two consecutive sites 2k — 1 
and 2k. At the same time the shock position lies between two consecutive sites 2k and 
2k + 1 in |/i2fc+i)- The authors in |8j have shown that the time evolution of the shock 
position in ([3]) governed by ([1]) is similar to that of a simple random walker. The shock 
position hops to the left and to the right with the probabilities 5i and 5 r respectively 
provided that: 

pi = , p 2 = 1 - ft 

Pi = ol , Pi = 1 (4) 
S r = f3 , Si = a. 

The steady-state of the system |P*) which obeys \P*) = T\P*) can be written as a 
superposition of the shocks defined in ([3]) as follows: 



^ 2n+l 
k=l 



\P*} = % E c*|A*> (5) 



c 2 fc = £ s (-^-) for k = -,n (6) 



provided that: 

• 5 Z ^2(n-fe)+l 
/(5;\2(n-fe) 

c 2 fc+i = (^-J for A; = 0, • • • , n (7) 

in which 5 S = 1 — 5 r — 5i. The partition function of the system Z (which is also the 
normalization factor) is also given by: 

In the steady-state and depending on whether 5i is larger or smaller than S r , the system 
undergoes a phase transition. The coexistence line occurs where 5i = 5 r . In both phases 
the density of particles on the lattice have exponential behaviors while on the coexistence 
line it is linear because of considering a superposition of shocks with unbiased random 
walk dynamics. 

In [6j the author has rigorously shown that the steady-state of the ASEP in a sublattice- 
parallel update can be obtained using a matrix product method recently reviewed in 
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[7]. According to this method we assign an operator to each state of a lattice site in 
the steady-state. The two operators D and E are associated with the presence of a 
particle and a vacancy at an even lattice site respectively. The two operators D and E 
are also associated with the presence of a single particle and a vacancy at an odd lattice 
site respectively. For a system of length 2n we associate a product of 2n operators to 
each configuration of the system consisting of particles and vacancies at different lattice 
sites. The unnormalized weight associated with any configuration can be obtained by 
considering a matrix element of this product. 

Using the matrix product method and in the same basis mentioned above one can rewrite 
the steady-state of the ASEP as follows: 



\n = ^((w\ 




|V» (9) 



in which the normalization factor is given by: 

Z = ((W\(E + D) n (E + D) n \V)). (10) 

It can be verified that the following two-dimensional matrix representation generates 
the same steady-state introduced in 



d= ? ,: u , e 





do 

fr(i-A0 






((W\ = (w u w 2 ), \V)) = 
given that: 

V2, (l-6 r )(p 2 - Pi) fS^ _ (1 -Si)(p 2 -pi) f . 

Vi (1 — d|) V (V W 2 (I -0 r ) 

Note that the two-dimensional matrix representation ( TTTT) also generates the same 
partition function (j5J) using (flQl) . In comparison to the matrix representation introduced 
in [8] we have adopted an upper triangular matrix representation. 
In what follows we will define a one-transit walk model consisting of 4n steps on a 
diagonally rotated square lattice and explain how the partition function of this system 
can be related to that of the ASEP studied above. We will then discuss how the physical 
quantities in both systems are related through a similarity transformation. 
Consider a random walker moving on a two-dimensional diagonally rotated square 
lattice. The random walker starts from the origin (0, 0) and takes 4n consecutive steps 
until it gets to the point (An, 0). The random walker moves in the north-east (NE) or 
the south-east (SE) direction. The random walker never takes two consecutive steps in 
the NE direction along the path; however, it can take two consecutive steps in the SE 
direction only once during the journey. In this case the random walker moves along a 
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Figure 1. Sketch of a simple one-transit path. 



path which crosses the horizontal axis only once. A typical one-transit path is given in 
Figure ((!]). Let us assume that the random walker moves along a weighted path. We 
assign a fugacity to the steps that the random walker takes as follows: 

• The fugacity z\ to each step in SE direction and the fugacity z 2 to each step in NE 
direction if the ending point is (2k, 0) for an odd k (k — 1, 3, • • • , 2n — 1). 

• The fugacity z\ to each step in SE direction and the fugacity z 2 to each step in NE 
direction if the ending point is (2k, 0) for an even k (k = 2, 4, • • • , 2n). 

• The fugacity 1 to all other steps. 

In order to find the weight of a given path one can simply multiply the fugacities of 
different steps of the path. The partition function of this system can be easily calculated 
and one finds: 

Z(z 1 , zi, z 2 , z 2 ) = (z 1 z 1 z 2 z 2 ) n Z(z 1 , zi, z 2 , z 2 ) (13) 

in which: 

Z(z 1 ,z 1 ,z 2 ,z 2 ) = ( jg ± - ( jg ± ^ (z 2 z 2 y\ (14) 

Z\Z\ Z2Z2 2121 2222 

The first term in ( IT3l) does not play any important role in the critical behavior of the 
system. The second term in f|T3|) . which is given explicitly in (I14p. can be reinterpreted 
as the partition function of a two-dimensional walk model on a diagonally rotated square 
lattice; however, one should assign the fugacities to the contact points (instead of the 
steps) in a different way. We assign a fugacity l/z\ to each contact point (2k, 0) for 
an even k and a fugacity 1/zi to each contact point (2k, 0) for an odd k if the contact 
point is above the horizontal axis. We also assign a fugacity l/z 2 to each contact point 
(2k, 0) for an even k and a fugacity l/z 2 to each contact point (2k, 0) for an odd k if 
the contact point is below the horizontal axis. In all cases we do not assign any fugacity 
to the contact points with the horizontal axis if they come from the first upward step 
or the first downward step. The partition function of this walk model given in (fT4"l) can 
now be written using a transfer matrix method as follows: 

Z(z 1 ,z 1 ,z 2 ,z 2 ) = (L\T n \R) (15) 

in which the transfer matrix T is a two-step transfer matrix defined as T = T Q T e . Note 
that T Q and T e are each related to two consecutive steps i.e. T Q (T e ) is the transfer matrix 
of two consecutive steps ending to the point (2k, 0) for an odd (even) k; therefore, one 
can rewrite T a and T e as: 

T = T° T e , T e = T° e T e e . (16) 
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Figure 2. The base vectors associated with the vertices and contacts. 



By assigning four base vectors to the vertices and contact points of a path one can find 
a matrix representation for the transfer matrices and the vectors introduced in (|15|) and 
([IS]) . The vectors \sk) and \s' k ) for k = 1, 2 are assigned to the vertices and the contact 
points of the path as shown in Figure ([2]). In an appropriate basis one can write: 

\si) = K> = ( I ) , h> = 14) = ( [ ) (17) 

with the following completeness property: 

E \*k)(8k\= E 14X41=2 (is) 

fc=l,2 fe=l,2 

in which X is a 2 x 2 identity matrix. Considering the definition of the weighted walk 
model the only nonzero matrix elements are: 



(si\T e \s' x ) = 


j_ 


(si|T e e K) = 


j_ 

z\ 


<*M«i> = 


1 




1 


<*M*a> = 


1 


(«i|I?|aa> = 


1 


(s2\T e \s' 2 ) = 


J_ 

Z2 


(^ 2 |T e e | S ' 2 ) = 


J_ 

Z2 


W 2 \T °\s 2 } = 


1 


, W2\T°\s 2 } = 


1 


(s[\R) = 1 




(s' 2 \R) = 1 




(LWi) = 1 




(L\s' 2 )=0 





All other matrix elements are zero. Now using f lT9|) one can easily construct the transfer 
matrices in the above mentioned basis: 

T °° = Te ° = (o i)' T ° e= (t °)' Te£= (o 1) (20) 

and the vectors (L\ and \R) are: 



(L\ = (l ) , \R) = 




Before going any farther let us investigate the asymptotic behavior of the partition 
function of the one-transit walk model given in ( I14p in the large-n limit. It can be 
easily seen that in the large-n limit the partition function shows two different behaviors 
depending on whether Z\Z\ > z 2 z 2 or z\%\ < z 2 z 2 which results in two different phases. 
This means that the phase diagram of the walk model depends only on z\%\ and z 2 z 2 
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and not the four independent fugacities. It turns out that for the case Z\Z\ > z 2 z 2 most 
of the contacts with the horizontal axis will be from below while in the case Z\Z\ < z 2 z 2 
they will be from above. On the coexistence line where z±zi = z 2 z 2 the number of 
contacts both from above and below the horizontal axis will vary linearly along the 
horizontal axis. 

At this point we show that how the partition function of the walk model fTHj) is related 
to that of the ASEP given in ([8]). It will also leads us to understand how the physical 
quantities in these systems are related. Let us assume that: 

Z\Z\ = 1, z 2 z 2 = (-^) 2 , z x z 2 = - p -. (22) 

Ol l — 0r — Oi 

Now the following similarity transformation connects the two above mentioned systems: 
((W\U=(L\, lf- 1 C 2 U = T, U- 1 \V)) = \R) (23) 

in which: 

It is also necessary to have the following relations: 

l — 5 r f5i\ _! l — 6i 



Wl = l w 2 = y—l _W ,vi = 7 - Y ,v 2 = u. (25) 

provided that V\, v 2 , W\ and w 2 satisfy (|T2|) which is achievable since u, d and d are 
free to be chosen. It is now easy to check that the partition functions of the two systems 
in different descriptions given by (jH]), ffTUj) and f JT4|) are equal. 

Let us now investigate how the physical quantities in the walk model are related to 
those in the ASEP. We start with the probability of finding a contact at the lattice site 
(2k, 0) for an even k from above P 2 j^ ove or below P 2 f ow the horizontal axis which is given 
by the following expressions: 

PS- = (LlT l^!^ m k = 2, 4, 6, ■ • • , 2„, (26) 

tit" = (UT \l^ R ) W * = 2, *, 8, • • • , ». (27) 

in which A 2 k = |si)2fc-i 2k( s i\ an d ^2fe = 1 52)2*1-1 2k(s 2 \ are the contact operators at 
the lattice site (2k, 0) from above and below the horizontal axis, respectively. Using the 
matrix representations (|20|) and (I2~T|) . it is easy to calculate the expressions (1261) and 
(1271) analytically. 

One can also calculate the probability of finding a contact at the lattice site (2k, 0) for 
an odd k from above P 2 l° ve or below P^ low the horizontal axis which is given by the 
following expressions: 

*ST = {L]T ^ T ^K) ^ m * = 1. 3. S. • • • , 2» - 1, (28) 
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#3f— = (L|r - T [^\I)~ ' ' fl> * = 1. s, 8, • ■ ■ , a» - 1. (29) 

The contact operators A 2 k and i?2fc have the same definitions mentioned above. 
Now we show that how the density of the particles in the ASEP is related to the 
probabilities of contacts given in ( )26]) -( l29]) . According to the matrix product formalism 
the density of the particles in an even site (r^) (k even) in the ASEP in the sublattice- 
parallel update is given by: 

Using ( 1231) . ( |26|) and ( 1271) and after some straightforward calculations it turns out that 
( 130]) can be written as: 

(r k ) = Pl P 2 t ve + P2P\t° w for k = 2, 4, 6, • • • , 2n. (31) 

On the other hand, in the same system and using the matrix product formalism the 
density of the particles in an odd site (r^) (k odd) is given by: 

w JW ( ?S» V " fat - w "- > '- 1 ' (32) 

Using ( 12"B"|) . (I2"gj) and (I2"§j) it is not difficult to show that (I3"2"j) can be written as: 

(r k ) = PxPT e + foffir for k = 1, 3, 5, • • • , 2n - 1. (33) 

In (5j E| the authors have already studied both the phase diagram and also the behaviors 
of ( )3T|) and ( 1331) in the large-n limit; therefore, we do not repeat their results here and 
the reader can directly refer to those references. On the other hand, a description of the 
phase diagram of the same system based on the traveling shocks has been done in [8]. 
In this paper we aimed to show that the steady-state of the ASEP with open boundaries 
in a sublattice-parallel update can be described using a two-dimensional walk model. 
We showed that the partition functions of these models are equal by properly defining 
the fugacities assigned to the steps (or the contact points) of the random walker. On 
the other hand, the probability of finding a particles at a given lattice site in the ASEP 
is connected to the probability that the Dyck path touches the horizontal axis at that 
point. This is valid for both even and odd lattice sites as obtained in ( l3Tj) and (133]) . 
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